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Abstract: Salt-affected soils are mostly found in irrigated areas within arid and semi-arid regions where 
the groundwater table is shallow. Soils of this type have become an increasingly severe problem because 
they threaten both the environment and the sustainable development of irrigated agriculture. A tool to 
estimate phreatic evaporation is therefore urgently required to minimize the salinization potential of 
salt-affected areas. In this context, phreatic evaporation at zero water table depth (Eo) is a key parameter 
for establishing a model for calculating phreatic evaporation. The aim of this study was to explore the law 
of phreatic evaporation and to develop structurally rational empirical models for calculating phreatic 
evaporation, based on Eo data of six types of soil (i.e., gravel, fine sand, sandy loam, light loam, medium 
loam, and heavy loam) observed using the non-weighing lysimeter and water surface evaporation (E601) 
data observed using a E601 evaporator of same evaporation area with a lysimeter-tube at the groundwater 
balance station of the Weigan River Management Office in Xinjiang Uygur Autonomous Region, China, 
during the non-freezing period (April to October) between 1990 and 1994. The relationship between Eo 
and Foo1 was analyzed, the relationship between the ratio of Eo to Foo and the mechanical compositions 
of different soils was presented, and the factors influencing Eo were discussed. The results of this study 
reveal that Eo is not equal to Foo. In fact, only values of the former for fine sand are close to those of the 
latter. Data also show that Eo values are related to soil texture as well as to potential atmospheric 
evaporation, the ratio of Eo to Foo and the silt-clay particle content (grain diameter less than 0.02 mm) is 
negatively exponentially correlated, and that soil thermal capacity plays a key role in phreatic evaporation 
at Eo. The results of this analysis therefore imply that the treatment of zero phreatic depth is an essential 
requirement when constructing groundwater balance stations to study the law of phreatic evaporation. 
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1 Introduction 


Phreatic evaporation from bare soil is a process that transfers water from phreatic zone to 
unsaturated zone before it passes into the air via soil evaporation. Understanding this process is of 


“Corresponding author: HU Shunjun (E-mail: hushunjun @ms.xjb.ac.cn) 

Received 2017-12-27; revised 2018-03-29; accepted 2018-05-14 

© Xinjiang Institute of Ecology and Geography, Chinese Academy of Sciences, Science Press and Springer-Verlag GmbH Germany, 
part of Springer Nature 2018 


http://jal.xjegi.com; www.springer.com/40333 


HU Shunjun et al.: Comparing phreatic evaporation at zero water table depth with water surface evaporation 


great importance to a range of research fields, including hydrology, hydrogeology, soil physics, 
and agriculture, especially in the context of evaluating water resource and calculating farmland 
drainage (Jin, 1982; Qu et al., 1983; Zhang, 1986, 1996, 2013; Huo and Wang, 1988; Lei et al., 
1988; Shen, 1992; Wang, 1993; Wang et al., 1994; Zhang and Shen, 1998; Ma et al., 2006; Su 
and Tao, 2009). In arid and semi-arid regions where areas with a shallow groundwater table, this 
process is a particular issue because considerable amounts of water are lost to the atmosphere and 
salts are accumulated within upper parts of the soil profile. In places where the groundwater table 
is shallow, evaporation can also lead to increasingly severe soil salinization problems which result 
in a surface salt crust that hampers agricultural activities (Zarei et al., 2002, 2009); such crusts 
also threaten the environment and sustainable agricultural development in these regions 
(Williams, 1987) and are therefore critical factors affecting the stability of oasis agro-ecosystems. 
Soil secondary salinization is particularly serious within the Xinjiang region of northwestern 
China because of the evaporation of groundwater with less rainfall, a larger evaporation rate 
overall, and a shallow water table (Shi et al., 2013). At the same time, however, direct evaporation 
measurements over large scales are problematic, costly, time consuming, and often impractical 
(Zarei et al., 2010). It is therefore essential to establish a tool to estimate phreatic evaporation 
(capillary rise) to minimize the salinization potential within salt-affected areas. In this context, 
Averianov (1956) proposed an empirical equation to estimate phreatic evaporation (Eq. 1). 


H n 
E =B- (1) 
where Eg denotes phreatic evaporation intensity (mm/d) when the water table depth (m) is H (i.e., 
from the ground surface to the groundwater table); Eois the phreatic evaporation intensity (mm/d) 
when the water table depth is zero (a value which generally approximates surface water 
evaporation); Hmax is the critical water table depth (m) at which evaporation ceases; and n is the 
phreatic evaporation index. 

Even though the relationship between Eg and H is limited in expression, it has nevertheless 
been widely applied because it is a simple expression and requires fewer parameters (Zhang, 1986; 
Zammouri, 2001; Yang et al., 2011; Chen et al., 2013). At the same time, although the intensity 
of phreatic evaporation when the water table reaches the ground surface is an important parameter 
within the Averianov equation, limited research has been carried out regarding Eg when H is zero. 
There are currently four main views regarding the relationship between phreatic evaporation and 
water table depth. The first view is that the phreatic evaporation at zero water table depth (Eo) is 
equal to, or approximate to, water surface evaporation (E601) (Huo and Wang, 1988; Tang et al., 
1989; Zhang, 1996; Guan, 2010; Xu, 2011; Shen and Huang, 2015). The second view is that 
phreatic evaporation at Eo is more than Eo (Shi and Liang, 1964; Liang, 1992; Tindall et al., 
1999; Fang, 2009; Han et al., 2009; Rui, 2014), while the third argument holds that phreatic 
evaporation at Eo is not equal to Feo: (Zheng, 1951; Don and William, 1958). The final view in 
this context suggests that phreatic evaporation is the same as E60: when the depth of the water 
table is less than decoupling depth or that of the evapotranspiration (ET) surface (Wang, 1993; 
Cook and Rassam, 2002; Shah et al., 2007; Ridolfi et al., 2008; Johnson et al., 2010). In an earlier 
study, Saxena et al. (1971) noted that water is lost readily from most soils at rates comparable to 
those from a free water surface when the water table is shallower than 60 cm, while Hillel (1998) 
considered Eo to be equivalent to potential evaporation. In a later study carried out in Tunisia, 
Zammouri (2001) noted that the Averianov equation provides more accurate estimates of Eg than 
either the Gardner or MODFLOW approaches, while Blight (2002) concluded that more 
evaporation may occur from the bare soil surface than from the water surface. In a series of 
studies carried out within the Tarim River Basin, Yu (1993), Mao et al. (1998), Shang et al. (1999), 
and Hu et al. (2009) established a series of empirical equations for estimating E, based on 
experimental data collected at the Qongkol, Shache, and Aksu water balance experimental 
stations in South Xinjiang, but in the absence of Eo measurements. Additional information 
regarding phreatic evaporation intensity when the water table depth is zero should therefore be 


JOURNAL OF ARID LAND 


obtained to enable further studies regarding whether, or not, the phreatic evaporation at Eo 
approximates F¢01. 

The objectives of this study are to: (1) identify whether phreatic evaporation for different soil 
textures is equal to E601 (as measured by the E601 evaporator) when water table depth is zero; (2) 
analyze the relationships among phreatic evaporation at Eo, E601, and the mechanical composition 
of different soils; and (3) provide a parameter to establish a calculation model for phreatic 
evaporation. 


2 Materials and methods 


2.1 Study site 


The experiments described in this study were conducted using a non-weighing lysimeter (Fig. 1) 
at the groundwater balance station (40°30’N, 82°33’E; 1013 m a.s.l.) operated by the Weigan 
River Management Office in Xinjiang Uygur Autonomous Region, China. The annual mean 
temperature is 10.7°C, mean annual precipitation is 53.2 mm, mean annual water surface 
evaporation is 2400 mm measured by an E®20 evaporator, and annual frost-free days are 216 d. 
The maximum frozen soil depth is 80 cm, and there are 2900 annual hours of sunshine. The 
Weigan River Management Office experimental station encompasses 120 hm? where the soil 
mostly consists of a grey sandy loam and a yellow mild clay; these clay layers occur 2 m below 
the land surface within the profile. 
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Fig. 1 Schematic illustration of the non-weighing lysimeter used in this study 


Drainage collector 


2.2 Experimental method 


To accurately encapsulate the range of soil textures seen within the Weigan River irrigation 
district, we filled a series of 54 lysimeters (cross-sectional area, 0.3 m?) with six types of soil 
(gravel, fine sand, sandy loam, light loam, medium loam, and heavy loam) at nine water table 
depths (0.00, 0.25, 0.50, 0.75, 1.00, 1.50, 2.50, 3.50 and 4.50 m) to observe phreatic evaporation 
(Fig. 2). The mechanical compositions of these six soil types are summarized in Table 1 (Lai et 
al., 2003). The lysimeters used for these experiments comprise a phreatic evaporation measuring 
cylinder, a Marriotte bottle, a drainage collector, and a connector (Fig. 1). Cylinders used to 
control water table depths less than 1.00 m were installed on both sides of the observation gallery, 
while those used to control water table depths more than 1.00 m were installed around the 
underground observation room (Fig. 2). A Marriotte bottle system was connected to the soil 
column in each case to control and record the internal water table and to measure the amount of 
water supplied to the soil column (Fig. 1). Gravity drainage was recorded by the collector, while 
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the water supplied to each lysimeter from the Marriotte bottle was measured three times a day at 
08:00, 14:00, and 20:00 (LST). Meteorological data were measured at a standard weather station 
within the experimental area; variables included air and ground temperature, humidity, hours of 
sunshine, wind speed at 2 m above the ground, rainfall, vapor pressure, and E601. 
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Fig. 2 Schematic illustration of the Weigan River groundwater evaporation station 


Table 1 Soil mechanical composition 


Weight percentage of soil particle content (%) 


Soil texture 


<0.02 mm 0.02-0.25 mm >0.25 mm 
Gravel 0.00 0.00 100.00 
Fine sand 6.48 37.07 56.45 
Sandy loam 32.36 67.11 0.53 
Light loam 40.28 54.48 5.24 
Medium loam 41.79 54.44 3.77 
Heavy loam 54.95 44.22 0.83 


3 Results 


3.1 Relationship between phreatic evaporation at zero water table depth (Eo) and water 
surface evaporation (E601) 


Data reveal the relationship between the phreatic evaporation intensity of the six soil types in the 
Weigan River irrigation district when water table depth is zero and E¢o: (Fig. 3). Results reveal 
that the phreatic evaporation intensity of all six soil types at Eo within this region is linearly 
related to Feo: (Fig. 3). This relationship is extremely significant as the correlation coefficient is 
greater than 0.86 and can be expressed as Equation 2. 

Eo= B x Eo s 2) 
where p denotes the slope of the correlation equation between phreatic evaporation intensity when 
water table depth is zero and E601 (mm/d), while E601 denotes the water surface evaporation 
measured by the E601 evaporator. 

Thus, if J equals 1, then Eo is equal to E601. Results further reveal the slope of the correlation 
equation between phreatic evaporation of the six soil types when water table depth is zero and 
Eoo1 within the Weigan River irrigation district (Table 2); as these data show that the £ value 
corresponding to the six soil types does not equal one, replacing Eo with E60: and establishing an 
empirical equation to calculate phreatic evaporation corresponding to different water table depths 
will generate a larger error. 
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Table 2 Correlation equation slopes between Eo and E601 


Soil texture Slope £ R R Roo Significance level 

Gravel 1.2778 0.9257 0.9621 0.2944 Extremely significant 
Fine sand 1.1026 0.8866 0.9416 0.3204 Extremely significant 
Sandy loam 0.7690 0.7538 0.8682 0.3395 Extremely significant 
Light loam 0.6606 0.7397 0.8601 0.2798 Extremely significant 
Medium loam 0.5586 0.7600 0.8718 0.3089 Extremely significant 
Heavy loam 0.4999 0.7669 0.8757 0.3612 Extremely significant 


Note: Eo, phreatic evaporation at zero water table depth; E01, water surface evaporation; J, slope; R’, determination coefficient; R, 


correlation coefficient; Roo1, critical correlation coefficient. 
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Fig. 3 Relationships between phreatic evaporation at zero water table depth (Eo) and water surface evaporation 
(E601) for six soil types 
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3.2 Relationships of Eo, E601 and soil mechanical composition 


The data presented in Figure 4 illustrate the relationship among Eo, E601, and the percentage of 
silt-clay particle (<0.02 mm) content (C) of a soil type. The ordinate in this figure is the ratio 
between Eo and E601, the slope of the correlation equation between phreatic evaporation intensity 
when water table depth is zero and E01. This was obtained by regressing Equation 2 (Table 2), 
and can be fitted with an exponential function (Eq. 3). 
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Fig.4 Relationship between Eo/E601 and silt-clay particle content 


4 Discussion and conclusions 


The Averianov equation shows that phreatic evaporation at Eo is equal to E601. The experimental 
results presented in this study reveal that Eo is greater than E601 in the case of gravel, but less than 
E601 in the case of loam, only approximate E601 in the case of fine sand within the Weigan River 
irrigated area. This analysis therefore suggests that additional research on phreatic evaporation at 
a water table depth of zero would be worthwhile. 

Phreatic evaporation at Eo refers to the evaporative process from a saturated soil surface (Jin, 
1982) as well as directly from phreatic water. As an oasis effect is also present when observing 
phreatic evaporation and Eo; at the same time, the evaporation area within phreatic evaporation 
measuring cylinders was held the same as in the E601 evaporator. Because environmental 
conditions were also the same throughout experiments, Equation 2 shows that applying the ratio 
between Eo and E601 weakens the oasis effect in this study. 

The experimental results presented here also show that phreatic evaporation at Eo decreases in 
concert with soil particle size. This corroborates the earlier findings of Don and Willian (1958), 
but conflicts with some other previous researches. It is generally recognized that when soils are 
very wet, the evaporation of water is governed by external atmospheric conditions rather than by 
their inherent properties; evaporation rate is therefore equivalent to potential rate, referred to as 
"potential evaporation", i.e., E601 ( Hanks and Ashcroft, 1984; Jury et al., 1991; Hanks, 1992; 
Hillel, 1998, 2004; Lal and Shukla, 2004). Phreatic evaporation at Eo is not only related to 
potential atmospheric evaporation, but also is influenced by evaporation surface roughness, soil 
heat capacity, soil thermal conductivity, soil color, and soil surface albedo. 

The data assembled in this analysis show that the particle sizes of the six soil types decrease in 
concert with evaporative surface roughness. Indeed, as an aerodynamically rough surface 
increases the turbulent mixing of air, Eo for gravel is greater than Feo. As a rough irregular 
surface also has a greater evaporation area (Tang et al., 1992; Wang, 2016), phreatic evaporation 
at Eo should increase (Wu and Zhang, 1986; Liang, 1992; Tindall et al., 1999; Cui, 2006; Lu and 
Pu, 2006; Fan et al., 2014; Shen and Huang, 2015). This theoretical result is inconsistent with the 
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experimental results of this study. 

Soil thermal capacity and thermal conductivity are closely related to mechanical composition. 
Thus, if soil thermal capacity is large, the temperature rise will be small when the surface is 
heated, water molecules will move slowly, and little water evaporation will occur. Previous 
research has also shown that soil heat capacity increases in concert with water content (De Vries, 
1963; Jury et al., 1991; Moldup et al., 2000; Wang, 2016). This means that as soil particle size 
decreases, saturated soil water content increases (Cosby et al., 1984) alongside thermal capacity 
while phreatic evaporation at Eo should decrease (Hu and Hou, 1988; Fang, 2009; Rui, 2014). 
This theoretical result agrees reasonably well with the experimental results of this study. 

The net radiation (Rn) received by a bare soil surface is partitioned into heating the soil and 
atmosphere, and evaporating water. Thus, under steady-state conditions (i.e., time-invariant and 
horizontally homogeneous), surface energy budget can be expressed as follows (Tindall et al., 


1999): 
Ri=S+LE+G, (4) 


where S denotes sensible-heat flux (W/m2); LE denotes latent-heat flux (W/m2); and G refers to 
soil heat flux (W/m?) (Tindall et al., 1999). 
The heat flow equation can be written as follows (Hanks, 1992): 
G=-K, ae K, fo iv 
Az Az 
where Kq denotes soil thermal conductivity (W/(m-°C)); To denotes soil surface temperature (°C); 
and T; is the soil temperature at the soil depth of Az (m). 

This energy balance relationship predicts that when soil thermal conductivity is large, more 
heat will be conducted from the soil surface into deep soil, while latent heat will be lower and the 
amount of evaporated water will be small. This means that as soil particle size decreases, 
saturated water content increases (Cosby et al., 1984) and the thermal conductivity of a saturated 
soil will fall (Chung, 1987; Jury et al., 1991; Wang, 2016; Wang et al., 2016). Thus, in theory, 
phreatic evaporation at Eo should increase, an outcome that is inconsistent with the experimental 
results of this study. We therefore conclude that soil thermal capacity exerts a significant 
influence on phreatic evaporation at Eo. 

The potential phreatic evaporation refers to the removal of water from a fully saturated soil 
surface when the depth of the water table is less than the decoupling depth, the depth of ET 
surface, or the thickness of the capillary saturated zone. However, as upward water flow must be 
high enough to match the external rate, small variations in the position of the water table around 
zero will not markedly influence the evaporation rate so long as these changes remain less than 
the decoupling depth, the depth of ET surface, or the thickness of the capillary saturated zone. 

When all soil particles are fully submerged as a result of controlling the water table depth, it is 
clear that measured phreatic evaporation will be equal to E6o1. This phenomenon may explain 
why some researchers have mistakenly noted that phreatic evaporation at Eo is equal to E601. 

Soil color affects surface albedo and the absorption of solar radiation by the soil (i.e., soil color 
influences soil evaporation by changing evaporation surface temperature). In other words, the 
darker the soil, the more the heat is absorbed by soil, the faster the soil surface temperature rises, 
the higher the evaporation surface temperature, and the greater the evaporation. Experimental 
data show that evaporation from a brown soil is 19% higher than that from a white soil, while 
evaporation from a black soil is 32% higher than that from a white one. Similarly, evaporation 
from a yellow soil will be 7% higher than that from a white soil (Shi, 1959; Fan et al., 2014; Shen 
and Huang, 2015). The influence of soil color and albedo on phreatic evaporation at Eo also 
require further study. 

The results of this study lead to the suggestion that phreatic evaporation at Eo should be taken 
into account when groundwater balance stations are constructed to study the phreatic evaporation. 
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